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ABSTRACT 

We present X-ray data on the distant cluster C10016+16 (z=0.5545) from ROSAT 
PSPC and HFJ observations and use them to study the physics of the intracluster 
medium (ICM) and the dynamical state of the cluster. The surface brightness dis- 
tribution is not only described by a spherically symmetric model but also by a two- 
dimensional (3- model fit. Subtracting an elliptical model cluster as defined by the best 
fit parameters of the two-dimensional model we find significant residuals, indicating 
an additional, extended X-ray source within the cluster. This source, likely to be a 
merging subcomponent of the cluster, coincides with a peak in the weak lensing mass 
map of Smail et. al. (1995). In the course of this analysis we present a new approach to 
quantify the significance of substructure in cluster X-ray images dominated by Poisson 
noise and smoothed with a Gauss filter. 

We determine the radial mass profile integrated out to a radius of 3Mpc and find for 
the total mass of the cluster a value of ~ 1.4 - 3.3 x 10 15 M Q and ~ 4.5 x 10 14 M for 
the gas mass, yielding a gas-to-total mass ratio of 14 - 32%. There is no significant 
radial dependence of the gas-to-total mass ratio in the cluster. 

Key words: galaxies: clusters: individual: C10016+16 - intergalactic medium - grav- 
itational lensing - X-rays: galaxies - cosmology: dark matter 



1 INTRODUCTION 

C10016+16 is with a redshift of z=0.5545 one of the best 
studied clusters of galaxies at higher redshifts. It has been 
target of observations in all wavelength. The cluster is very 
massive and, comparing it to nearby clusters, it certainly is 
most similar to the Coma cluster. It seems to be embedded 
in a large-scale-structure density enhancement at a redshift 
of approximately z=0.55 (Koo, 1981). This idea was recently 
strengthened by Hughes, Birkinshaw & Huchra (1995), who 
found a poor cluster in X-rays at a redshift of 0.5506 at 
about 8 arcmin distance from C10016+16. 
It is an exceptional cluster, due to its high X-ray luminosity 
(L a; (2-10keV)=2. 62 xlO 45 erg/sec (Tsuru et al. 1996) at such 
a high redshift. Also the optical appearance of this cluster is 
different to other equally distant clusters. For a long time it 
was known to be a counter example for the Butcher-Oemler 
effect (Butcher& Oemler 1978) showing only a red, old pop- 
ulation of galaxies. 

Belloni& Roser (1996) (hereafter BR) recently found that 
the cluster has a high fraction of E+A galaxies, being re- 
sponsible for the relatively high fraction of red light in the 
cluster's galaxies. Wirth, Koo & Kron (1994) studied 24 
cluster members with the HST, and found that the E+A 
galaxies seem to be more disklike than normal red galaxies. 
C10016+16 is an ideal cluster for the observation of the 



Sunyaev-Zel'dovich effect, as it is very X-ray luminous and 
at a high redshift. The cluster was among the first three 
objects for which a detection of the effect has been claimed 
(Birkinshaw, Gull & Hardebeck 1984; Uson 1986; see also 
Rephaeli 1995). 

In this paper we study the X-ray properties of C10016+16 us- 
ing ROSAT/PSPC and HRI data. We give constraints on the 
radial total mass profile of the cluster and discuss its mor- 
phology and its dynamical state. Most probably C10016+16 
is in the process of a merger. The smaller infalling compo- 
nent carries only a small fraction of the mass, and therefore 
the merger has probably a small effect on the overall dy- 
namical equilibrium of the cluster. 

Recently, a weak gravitational shear signal was detected in 
the cluster by Smail et al. (1995) (hereafter SEFE). They 
used this measurement to infer the gravitational mass of 
C10016+16. We compare our results on the cluster mass de- 
rived from the X-ray data with their weak gravitational lens- 
ing results. 

Throughout the paper we use a Hubble constant of 
H =50km /sec /Mpc. 



2 THE OBSERVATIONS 
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Figure 1. Contour Plot of the ROSAT HRI observation of 
C10016+16. The image is smoothed with a Gauss filter with a a 
of 6 arcsec. The lowest contour level is 1.3xl0 — 6 cts/sec/arcsec 2 . 
The spacing is linear with steps of 5.22xl0~ 7 cts/sec/arcsec 2 . 

2.1 ROSAT HRI data 

C10016+16 was observed by the ROSAT/HRI for 76593 sec 
(for a description of ROSAT see Triimper 1983, 1992). The 
observation was performed in two different periods, in Jan- 
uary and June/July 1995 with exposure times of 5,917 sec 
and 70,676 sec, respectively. To correct for a possible offset 
between the two pointings, we use the QSO 0015+162 in 
the North of the cluster. We measure an offset of 6" in right 
ascention and 4" in declination between both positions of 
the QSO, and correct for it. 

Fig.0. shows the HRI count rate image of C10016+16. The 
image covers the total energy range of the ROSAT telescope, 
0.1-2.4keV. 

The source in the North-East of the centre is not very sig- 
nificant, and does not have an optical counterpart (see also 
Fig.|.). The QSO 0015+162 lies outside of this plot but can 
be seen in Fig.nl 



2.2 The PSPC data 

C10016+16 was observed for 43,157 sec with the ROSAT 
PSPC (for more detailed information see Hughes et al. 
(1995)). To obtain the best signal-to- noise-ratio for the de- 
tection of the cluster we select the photons in the 0.5 - 2 keV 
energy band (channels 52 to 201). The resulting count rate 
image for this energy range is shown in Fig|| The image is 
vignetting corrected. The PSPC data show a much more reg- 
ular appearance, partly coming from the larger point spread 
function (PSF) of the instrument in comparison with the 
HRI, and partly from the better statistics. The PSPC has 
a much lower background and a higher sensitivity (about a 
factor of three) compared to the HRI. 

The most prominent point source in the image is QSO 
0015+162, lying to the North of the cluster, approximately 
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Figure 2. Contour plot of the exposure corrected PSPC data 
of C10016+16. The image is smoothed with a variable Gauss 
filter. The largest a is 30 arcsec. The lowest contour line is 
1 X 10 _7 cts/sec/arcsec 2 The spacing of the contour levels is 35% 
of the value of the lower one. 



3 arcmin from the centre. In the South, approximately 8 ar- 
cmin from the centre, one can see a clearly extended source 
which is another galaxy cluster at a redshift of z=0.5506 re- 
cently discovered by Hughes et al. (1995) in the same PSPC 
data. 



3 THE SPATIAL DATA ANALYSIS 

3.1 The spherical symmetric fitting 

For the data analysis we use EXSAS, the software system 
provided by MPE (Zimmermann et al. 1994). To study the 
global physical properties of the cluster we primarily use 
the PSPC data, because they trace the cluster X-ray emis- 
sion to larger radii due to the higher signal to noise and 
better photon statistics. But we also present the results of 
analyzing the HRI data. With the PSPC data we have also 
checked for a possible hardness ratio variation of the X-ray 
emission across the PSPC cluster image and found no vari- 
ation, implying that the soft energy band does not provide 
any particular extra information on the physical parameters 
of the cluster. 

To obtain an approximate description of the X-ray surface 
brightness distribution, which allows to subsequently de- 
rive the main physical parameters of the cluster analytically, 
we use the so-called isothermal /3-model (Cavaliere&Fusco- 
Femiano, 1976, 1981; Sarazin&Bahcall, 1977; Gorenstein et 
al., 1978; Jones&Forman, 1984), which describes the sur- 
face brightness Sir) of a galaxy cluster assuming spherical 
symmetry 

/ 2 s -3,3 + 1/2 

S(r) = So f 1 + ^ J + B (1) 
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where So is the central intensity, r the radius, r c the core 
radius, (3 a slope parameter of the radial surface brightness 
distribution, and B the background surface brightness. The 
X 2 -fitting results for the PSPC data are: (3 = 0.80j 
r c = 372±gkpc= 50.5± 8 7 j" (for a summary see Tab.j 
surface brightness profile and the fit are shown in Fig 
results for the HRI data are: r c = 283±^° kpc, (3 = O.i 
(see also Tab|j). Serendipitous sources are extracted from 
the fitting. The best fit values of both data sets for the core 
radius and the /3 are quite different, but the error bars have a 
large overlap (see Tab. [2] and Tab.g). The discrepancy arises 
from the fact, that the HRI data do not trace the cluster 
out to such a large radius as the PSPC data, due to higher 
background and lower sensitivity. 

Fitting this spherically symmetric profile to the image of 
C10016+16 has the disadvantage that it does not perfectly 
describe the slightly elliptical cluster. But it allows a sim- 
ple analytic deprojection of the surface brightness profile of 
the cluster yielding the radial gas density profile. We will 
show later that the spherically symmetric approximation is 
a sufficiently good description for the determination of the 
gas mass and total mass profiles with errors that are smaller 
than relevant uncertainties like those of the temperature dis- 
tribution of the intracluster medium (ICM). It can actually 
be seen also from the analysis of simulated clusters that 
a spherically symmetric /3 model approximation generally 
leads to a good description of the shape of realistic clusters 
with errors of less than 20% (e.g. Schindler 1995; Evrard, 
Metzler & Navarro 1996). The analytic deprojection of the 
profile of equ. (1) leads to the density profile: 



p(r) 



Po 



-3/3/2 



(2) 



For the central electron density we obtain 0.65x10 2 cm 3 
from the PSPC data, and 0.77x 10~ 2 cm~ 3 from the HRI 
data. 

For the fitting we cut out regions of serendipitous sources, 
which clearly do not belong to the emission of the ICM it- 
self. 

The different fit parameters are correlated. For example the 
j3 and the core radius are dependent on each other. But 
(3 is also correlated to the background. A too low result 
in the background leads to a decrease of f3, as the profile 
has to become shallower to overcome the difference of the 
wrong background result to the real one. Therefore we have 
carefully checked these sources of uncertainty and found in 
particular that our result for the background in the fit is in 
perfect agreement with the background determined indepen- 
dently in large regions outside the cluster with all detectable 
individual sources removed. 



3.2 The two-dimensional fit 

To account better for the slightly elliptical shape of the clus- 
ter we have also performed a two-dimensional fit to the clus- 
ter using a modified (3 model that allows for two different 
core radii along the two principal axes of the cluster image 
ellipse. The surface brightness profile in this model is de- 




Figure 3. Surface brightness profile of C10016+16 from the PSPC 
data. The crosses mark the errors and the full line is the best fit 
with the parameters shown in Tab.^. 



M 

Table 1. Description of the fit parameters of the two-dimensional 
case. Not mentioned parameters have the same definition as in the 
Id fit. 



parameter 



definition 



I'D 

2/0 

n 

a 
ai 

0,2 



central position in x direction 
central position in y direction 
coordinate in x direction 
coordinate in y direction 

position angle 
major axis of core radius 
minor axis of core radius 



scribed by the equation: 



S(r 1 ,r 2 )=S (l + F 1 +F 2 ) 



-3/3+1/2 



+ B 



(3) 



with 



Fi 



F 2 = 



(cos(a)(n - xo) - sin(a)(r 2 - yo)) 2 



(sin(a)(n - x ) + cos(a)(r 2 - yo)) 2 



The fit is applied to the two-dimensional pixel data of the 
image. The newly introduced parameters are defined in 
TabQ 

The results are shown in Tab.^ and Tab.^. This fit has 
the advantage of taking into account, that a relaxed cluster 
does not necessarily need to be completely spherical sym- 
metric, but can show ellipticity. A further advantage is that 
the centre position of the cluster is derived in the fitting 
procedure itself and does not have to be predetermined like 
it is necessary in the one-dimensional analysis. There the 
centre position is usually determined from the maximum in 
the surface brightness profile (with some dependence on the 
smoothing used before determining the centre). 
The application of such a two-dimensional fit can also be 
found in Bardelli et al. (1996) in which they apply this model 
to Abell 3558 the central cluster of the Shapley concentra- 
tion. 
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However, the two-dimensional fit exhibits more problems 
than the one- dimensional one, firstly because of poorer 
statistics in each bin, as one has to use pixels instead of 
concentric rings. This is particularly severe at larger radii 
because the surface brightness is falling off steeply, which is 
partly compensated by the increasing surface of the rings in 
the one- dimensional analysis but not in the two-dimensional 
case. This is very crucial in the regions outside the clus- 
ter emission, where the only observed emission comes from 
the background. Secondly, the number of fit parameters in- 
creases from four to eight fit parameters. 
For our two-dimensional fit we use the PSPC data with a 
binsize of 5" x 5". We again cut out serendipitous sources in 
the field of view. Since the number of photons per pixel in the 
background area is too small to apply Gaussian statistics, 
which is assumed by x 2 -fitting, we apply a small Gauss-filter 
to the image. The a of the Gauss-filter is 5", and therefore 
much smaller than the FWHM of the PSF (~ 30") of this 
data set, so that we do not risk to loose any information. 
If we do not apply a Gauss-filter, we underestimate the back- 
ground. This is due to the fact that for the low photon 
statistics in the background pixels the mean of the Pois- 
sonian distribution is larger than the most abundant value 
and the mean adopted by assuming a Gaussian distribution. 
This adopted Gaussian distribution is used with x 2 -fitting. 
Because of the above mentioned correlation effect an under- 
estimation of the background leads than to a decrease in the 
fitted slope parameter j3. The application of a small Gauss 
filter improves the photons statistics and removes this prob- 
lem. Of course, the errors obtained for the fitting parameters 
in this procedure are no longer precisely defined. However, it 
is a very good approach to get the best fit value. We test this 
with images on model clusters, on which we added Poisson 
noise. These model clusters have about the same parame- 
ters as C100I6+1 6. We apply a Gauss-filter, of the same size 
as for the real cluster image, and run the fitting routine on 
them. The parameters obtained from the fitting are in very 
good agreement with the intrinsic model parameters, so that 
we can be sure, that the Gauss-filter does not obscure our 
fitting result. 

The one-dimensional fit-parameters are in very good agree- 
ment with the two-dimensional fit- values (see Tab.^ and 
Fig.^), so that this can also be used as a proof for relia- 
bility of the Gauss filtering. The best fit value for the core 
radius in the one-dimensional fit is 372 kpc, lying almost ex- 
actly on the geometrical mean value of the 2d fit. Also the 
of the two fits agree very well, with 0.80 best fit value for 
the Id fit, and the 2d fit being at 0.81. 

Trying to overcome the Gauss filter with a larger binsize of 
the image is not feasible, due to the low background. The 
bins would be too large. The best fit value for from the 
two-dimensional fit are slightly depending on the binning. 
However, applying different binning and different Gauss fil- 
ters (not too large of course) yields values which are all in 
excellent agreement with the Id fit and its errors. 
Another approach, using a maximum likelihood analysis, 
as for example applied by Birkinshaw, Hughes & Arnaud 
(1991), has the disadvantage of predefining the fit parame- 
ters. The parameters are not fitted, but only tested. 
To check the existence of substructure, in addition to the 
overall ellipticity of the cluster, we produce a synthetic clus- 
ter image from the fit parameters of the two-dimensional 



model and subtract it from the real image. There is signifi- 
cant structure in the residual image that can indicate either 
point sources in the field or subcomponents of the cluster 
itself. The residual images are shown in Fig.^ for the PSPC 
data and in FigJH] for the HRI data. To obtain Fig|B| the 
HRI image was treated in the same way as the PSPC image 
with a two-dimensional fitting approach. The used parame- 
ters are shown in Tab.^. Here the pixel size is 4" x 4" and 
also a — 4" for the Gauss filter. The residual images shown 
in these figures are displayed in the form of a significance 
plot. How the significance of the pixels of the residual images 
are calculated is explained in more detail in the appendix. 
Briefly the initial problem is how to calculate the signifi- 
cance of a source after having applied a Gauss filter. We 
solve this problem by using error propagation and Poisson 
noise. Basically two different Gauss filtered images are di- 
vided by each other. The one being the "error image" , is not 
only Gauss filtered, but also the square root of each pixel 
has to be calculated together with a general normalization. 
The "source image" can be calculated in two ways. In the 
first case one subtracts a background model. Then the rou- 
tine gives a above background. The other possibility is to 
subtract a model cluster, and to calculate the significance 
(a) of these residuals. A more general approach for the test- 
ing of significances with filter techniques (not only Gauss 
filters) is described by Rosati (1995). 



4 THE SPECTRAL ANALYSIS 

Due to the spectral resolution of the PSPC it is possible 
to determine the temperature of the ICM. However, as this 
cluster has a very high redshift and only a relatively small 
number of source photons were detected, the spectral reso- 
lution is of course limited. 

We determine the temperature in different radii around the 
centre of C10016+16 with different backgrounds and differ- 
ent energy bands with different values for metallicity. The 
results show a large scatter, which is due to the relatively 
low photon statistics. For example in a radius of 4 arcmin 
around the centre, the total sum of photons is about 6,000. 
Nevertheless, the results which are in the range 6 to 10 keV 
are in quite good agreement with the ASCA data of Tsuru 
et al. (1996) who obtained a global temperature of 8. 22^47 
keV. Only the results for the hydrogen column density are 
different. We obtain results of 4 — 6 x 10 20 cm -2 , which are in 
good agreement with the measurements of the 21cm line of 
4 x 10 20 cm" 2 (Dickey & Lockman 1990), while Tsuru et al. 
find 12.9 x 10 20 cm -2 . This difference might be explained by 
the fact, that ASCA does not provide a very good spectral 
sensitivity below 0.5keV, which is the energy range relevant 
for the determination of this degree of absorption. 
In different rings (which are partly overlapping) we do not 
detect any significant temperature gradient, which leads us 
to the assumption that the cluster is probably more or less 
isothermal within the given the large error bars. This is in 
very good agreement with the ASCA data. Only in the very 
centre we see a drop of the lower boundary of the temper- 
ature (in a radius of 1 arcmin) to 4 keV together with an 
increase of the hydrogen column density. This decrease is, 
however, not significant and most likely an artifact of the 
temperature fitting, as the cluster has not yet developed a 
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Figure 4. Significance of the residuals after subtracting the el- 
liptical isothermal /3-model from the original hard PSPC image. 
The parameters for the subtracted model is described in Tab. El 
The stepsize of the contours is 1 a. The dashed line is the level, 
the full lines are positive the positive it's, the dotted ones are 
the negative ones. The numbers are the objects of the optical 
identification list of BR. The Gauss filter has a a of 25 arcsec. 




offset in RA in [arcmin] 

Figure 5. Significance of the residuals after subtracting the el- 
liptical isothermal /3-model from the original HRI image. The pa- 
rameters are presented in Tab|j. The stepsize is 1 <r. The dashed 
line is the level, the full lines are positive the positive <r's, the 
dotted ones are the negative ones. The numbers are the objects 
of the optical identification list of BR. The Gauss filter has a a 
of 20 arcsec. 



Cooling Flow (see also 5.2) for which such a temperature 
decrease would be expected. Therefore we neglect this lower 
boundary in the centre for the mass determination. Includ- 
ing it does not change the overall mass result significantly 
(the mean result changes by less than 5%). 
In general it is interesting, that the temperatures determined 
from the PSPC data are only very weakly depending on the 
assumed values for the metallicity. 

5 THE MORPHOLOGY OF CL0016+16 
5.1 The X-ray data 

C10016+16 shows a position angle of about 50° (measured 
from North to East) and an eccentricity of about 0.18. There 
are clearly additional sources superposed on the elliptical 
cluster, both in the HRI and the PSPC observation. 
In the West of the cluster, about 2 arcmin from the centre 
(FigjiJ), we find a surface brightness excess with 4 a sig- 
nificance in the PSPC image. This source at a position of 
RA=00H18M27S DEC=16D25M55S (J2000) is most likely 
extended, as explained in the discussion. 
The residual image of the HRI observation also shows a sur- 
face brightness excess to the West of the centre. At the exact 
location of the PSPC maximum of FigJ^ we find, however a 
gap in the HRI surface brightness excess (FigJEJ). There are 
several effects that can produce such a result. First of all, 
the HRI data have a lower photon statistics than the PSPC 
data as explained above and thus the gap may be produced 
by a larger statistical fluctuation (the surface brightness ex- 
cess is just produced by 55 photons). The second point is, 
that only the hard band photons were used for the analy- 



sis of the PSPC data. We cannot apply the same selection 
to the HRI, as this instrument does not provide any energy 
discrimination. But as we found no significant feature in a 
hardness ratio map produced from the PSPC data we be- 
lieve that the difference in the shape of the western emission 
excess is due to photon statistical effects. 
Therefore we conclude, that there is an extended source in 
the West in the line-of-sight of C10016+16. This extended 
source might be a subgroup falling into this cluster. How- 
ever, this subgroup, if existing is very small compared to the 
total cluster C10016+16. 

Another source visible in the significance plot of the PSPC 
image, but only having a 2a significance, is the source in the 
South. This source structure is most likely to be extended 
too. Galaxy number 43, defined by BR, having a redshift of 
z=0.38, indicated in the figure, coincides with the maximum 
of this source, so that it might be, that there is emission from 
a small foreground group, as already suggested by Ellis et 
al. (1985). However, optical source number 36 also in this 
enhancement is a star and might contribute partly to the 
emission. 

5.2 Cooling time 

The central cooling time, calculated for a gas temperature 
of 8.4keV (Tsuru et al, 1996; Yamashita, 1995) and for the 
density profile obtained from the isothermal f3 model fit is 
larger than 10 10 years from both the PSPC and the HRI 
data. Therefore the cluster cannot have yet formed a steady 
cooling flow (Fabian, Nulsen & Canizares 1984). 
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5.3 Comparison with optical data 

C10016+16 is a well observed cluster in the optical. There 
has been an extensive survey by Dressier & Gunn (1992) 
(hereafter DG) and BR. Figg and Figj| show the overlays 
of the optical data of the sample of BR. 
We did not apply a bore-side correction in the PSPC image, 
but found an inaccuracy of the data of 2 to 5 arcsec for the 
pointing position. 

The most striking feature is, that the alignment of the cen- 
tral galaxies is coinciding with the PA we determined from 
X-rays (as already mentioned by SEFE) . However, for the 4 
a excess emission in the West of the centre, which is proba- 
bly a substructure feature of the cluster, there is not enough 
overlap with the data of BR to search for a possible galaxy 
density enhancement which could be related to the surface 
brightness excess. DG detected objects in that region, but 
they did not obtain redshifts for them. The objects seem to 
be extended, so that it is likely that they are galaxies. 



6 MASS ANALYSIS 

6.1 The one-dimensional approach 

X-ray observations of the ICM in clusters of galaxies can be 
used for the determination of the total and gas mass. The 
formula for the cluster mass profile obtained from combining 
the isothermal j3 formalism (Id) equ(2) and the hydrostatic 
equation is: 



M(< r) = 



kr 2 



IdT 
\~d~r~ 



3/3rT 
r 2 + ri 



(4) 



Gfim 

The results are shown in Fig|j. For the mass determina- 
tion we use the PSPC data because of the better statistics, 
even though the spatial resolution is a factor about 7 worse 
in comparison to the HRI data. Our total mass estimate 
for C10016+16 at a radius of 3Mpc is about 2.3xl0 15 M Q . 
This is a typical result for a rich cluster, like for example the 
Coma cluster of galaxies (Briel, Henry & Bohringer 1992). 
For the mass determination we use the method by Neu- 
mann & Bohringer (1995). This method uses a Monte-Carlo 
approach for determining the cluster's radial mass profile. 
The stepwidth for the Monte-Carlo method is 300kpc with 
a maximal change of the temperature of 1 keV per step. 
For the temperature range we take 6 to 10 keV (rounded 
errors of Tsuru et al. 1996 and our results from the PSPC 
data fitting). We assume that the cluster gas is not neces- 
sarily isothermal, but only that the gas temperature lies in 
the given temperature range. This assumption is justified 
since we find no significant temperature variations in differ- 
ent concentric rings. The gas mass of the cluster integrated 
out to a radius of 3 Mpc is roug hly 4.5 xlO 14 M . Thus we 
obtain a gas to total mass ratio of about 20%. This is in the 
typical range of other clusters. 

6.2 The two-dimensional approach 

With our results for the two-dimensional fit, it is also possi- 
ble to determine the influence of ellipticity on the results for 
the cluster mass. For comparison we determine the masses 
with the major and minor-axis results for the core radius 



Table 2. Obtained parameters for C10016+16 using the PSPC 
data. The physical values are all calculated with Ho = 50 
km/sec/Mpc. The errors for r c and (3 from the Id fit are 3cr 
errors. 



parameter 



value 



/3 Id fit 
r c Id fit 
n e o Id fit 
Mtot < 3 Mpc 
M gas < 3 Mpc 



+0.09 
0.07 

+8.4// 
52 "-^ — ""'"-7.1 

0.65 xl0~ 2 cm- 3 



0.80Z 
372±62 kpc _ f)(| - 



2.321° ;g* x 10 15 M 
4.5 x 10 14 M 



(3 2d fit 

r c for major axis 
r c for minor axis 
PA (North over East) 
e (eccentricity) 



0.81 
410 kpc = 56" 
337 kpc = 46" 

50° 

0.18 



of the two-dimensional fit. We again use the Monte-Carlo 
approach with the same conditions as the one-dimensional 
analysis. As one can see, the difference between both results 
becomes less with increasing radius. This result is natural, 
as only the coreradius is different, which influences the mass 
result mostly near the centre. 



6.3 Comparison of the different mass estimates 

Comparing the different results for the total mass of the Id 
approach with the 2d approach leads to similar results, as 
can be seen in Fig.^. This also implies that it is a good ap- 
proach even to take only a sector out of a certain cluster 
surface brightness distribution when the cluster is ellipti- 
cal. This approach has been normally undertaken, if clus- 
ters show clear substructure, to exclude those regions. Only 
within the innermost 1 Mpc there is some divergence of the 
mass depending on the core radius. However averaging the 
2d distribution to a modified Id distribution leads to a core 
radius (by taking the geometrical mean of the major and 
minor axis) of less than 1% difference to the coreradius orig- 
inally derived by the Id fit. This is a proof, that the Id fit 
averages correctly over the eccentricity of a cluster. Also the 
similarity of the /3's of both approaches supports this. 
A comparison of the results for the gas mass shows that this 
case is slightly more critical. Because of the dependence of 
the central density on /3 and the coreradius, one yields dif- 
ferent results for the approach with the values for the major 
and minor axis. Therefore an analysis, that takes into ac- 
count only sectors of the X-ray distribution must be under- 
taken cautiously, not to overestimate or underestimate the 
mean coreradius. Taking the extreme values for the major 
and minor axis, we get a difference to the averaged value 
almost up to 20%. 

However, averaging over the entire cluster with a Id fit is a 
good approach as the averaged 2d central electron density 
is of the order of 1% different to the one from the original 
Id approach. Therefore the differences of the averaged Id 
approach to the original Id fit are in the order of less than 
1%, and only diverge to less than 1.5% at a radius larger 
than 2.5 Mpc. 
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Figure 6. Integrated mass profile of CL00f6+I6. The full lines 
indicate the mass profile from the Id fit. It shows the mean value 
and the upper and lower lines are the mean ±2<r. The small dots 
are the results for the 2d fit with the minor axis as core radius. 
The circles are the same for the major axis as core radius. Again 
middle values are the mean, the others ±2<r. The lower dotted line 
is the gas mass profile from the Id fit, the dashed line is the gas 
mass from the 2d fit, using the major axis value as core radius, 
the dotted dashed line respectively for the minor axis result used 
as core radius. 



Table 3. Obtained parameters for C10016+16 using the HRI 
data. The errors for r c and (3 from the Id fit are 3cr errors. 



parameter 



value 



P Id fit 
r c Id fit 

7l e Id fit 



0.68 



283 



+0.26 
-0.1 

kpc = 38.5 
0.76 xl0" 2 cm- 



-140 
-81 



+ 19/' 
11 



P 2d fit 

r c for major axis 
r c for minor axis 
PA (North over East) 
e (eccentricity) 



7 DISCUSSION 



0.69 

318 kpc EE 43" 
261 kpc = 35" 

50° 

0.18 



7.1 Possible influence of the QSO on the X-ray 
temperature 

Tsuru et al. (1996) found with ASCA a temperature of about 
8.22keV and a very low metallicity, m<0.167 for the ICM. 
However, these values can in principal be obscured as there is 
another X-ray emitting source, QSO 0015+162 in the vicin- 
ity of the cluster (see also Fig.|2j) . The distance to the cluster 
centre of 3 arcmin. This is in the order of the half power di- 
ameter of the point spread function of the ASCA plus CIS 
combination of 3.2 arcmin (Ikebe, 1996). Tsuru et al. 1996 
used the CIS data. It is very likely, that this QSO, which 
has a redshift of z=0.553, close to the cluster's redshift, in- 
fluences the spectrum and therefore the spectroscopic esti- 
mates might be wrong. A rough estimate of the influence 



can be made by comparing the countrate of the QSO with 
the countrate of the cluster itself measured by ROSAT. For 
the PSPC we get a countrate of about 10% for the QSO in 
comparison to the cluster. Fitting a power law to the QSO 
spectrum of the PSPC yields a power law index in the range 
of -3 to -2.3. Fitting a power law to the cluster itself yields a 
power law index of -1.6 to -1.5, both times fixing the hydro- 
gen column density to 5xl0 20 cm -2 . Thus the spectrum of 
the QSO is steeper implying that it has less influence in the 
higher energy band of ASCA. For the temperature measure- 
ment, the QSO would, if it influences the cluster spectrum, 
lower the fitted temperature, as fitting a Raymond-Smith 
spectrum to the QSO data yields a best fit parameter of 
about 2 keV. 

Another possibility in principle is to measure the influence of 
the QSO to the cluster spectrum by comparing luminosities 
measured with ASCA and the PSPC. Tsuru et al. (1996) ob- 
tain a value for L 2 ,(2-10keV) = 2.62 x 10 45 erg/sec. With the 
PSPC we get L x (2 ~ lOkeV) = 2.3 - 4.4 x 10 45 erg/sec. This 
shows that the errors are too large to give a definite result 
on the influence of the QSO. For the luminosity determina- 
tion we took into account variations of the hydrogen column 
density, the metallicity and the temperature (6-10keV). 
Taking this all together, it is not likely, that the QSO con- 
tributes sufficiently to the spectrum of the cluster, to explain 
the very low result for the metallicity. A 10% contribution 
to the continuum of the ICM spectrum probably does not 
lower the metallicity by a factor of two to three, which would 
be necessary, to make 00016+16 a 'normal' distant cluster, 
with a metallicity of about m=0.3. 

The low value for the metallicity is rather surprising, as 
C10016+16 is the only cluster showing such a result. This 
does not fit in the normal scenario of metal enrichment, es- 
pecially as the cluster exhibits many E+A galaxies (BR), 
which have already burned all O and B stars. These stars 
would have been in principle able to enrich the ICM by their 
supernovae. However, this does not seem to be the case. 



7.2 The dynamical state of C10016+16 

The X-ray image provides strong evidence for the existence 
of substructure in the cluster which is possible due to the 
merging of a galaxy group with the main cluster. 
In this case, there are basically two possibilities to explain 
the observation, as suggested by simulations (Schindler & 
Miiller 1993). The first one is, that the subgroup is on its 
first infall to the cluster, the second one is, that the sub- 
group already penetrated the cluster centre once. 
The second scenario leads to a much stronger distorted ap- 
pearance of the subgroup than the first one. As the subgroup 
can be seen as a more or less compact feature, the first case 
seems to be far more likely. Another very interesting feature 
is the chain of galaxies, which coincides very well with the 
PA of the gas. It is not clear whether this chain comes from 
a recent merger, and if, whether the measured substructure 
in the West has something to do with it. The merger direc- 
tion of the substructure roughly agrees with the PA of the 
gas. This is a common feature for mergers as predicted by 
simulations (Van Haarlem & Van de Weygaert, 1993). How- 
ever, it might also be, that the chain of galaxies are only 
the remnants of a former merger, proving that the cluster 
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formed out of filaments. 



7.2.1 Mass estimate of the possible subgroup 

The subgroup in the West shows up with about 55 photons 
in the hard ROSAT/PSPC band of 0.52-2.01 keV. 
Due to the limited number of photons, it is very difficult to 
give estimates on the physical quantities of this source. How- 
ever, assuming this source is at the cluster's redshift, and on 
its first infall to the cluster, so that it still contains most of 
its own X-ray emitting ICM, we obtain a X-ray luminosity 
of about L x ~ 10 43 erg/sec. Comparing this luminosity to 
groups of galaxies, we find, that this value is on the upper 
limit for Hickson groups (Ponman et al. 1996) . Poor clusters 
of galaxies, having more galaxies than Hickson groups are ly- 
ing in the regime of about 10 44 erg/sec, depending if they are 
hosts for Cooling Flows or not. Therefore we can classify our 
source as a group of galaxies in the range between Hickson 
groups and poor clusters of galaxies. The mass range, which 
one observes for these objects lies in the order of 10 13 Mq to 
10 14 M Q (Ponman et al. 1996; Mulchaey et al. 1996). This 
is the most likely mass range for the subgroup. Comparing 
this result with the total cluster mass of about 2xl0 15 M , 
we find a contribution of the subgroup to the total mass of 
about 1% to 20%. 



7.2.2 Contamination by stars 

It must be noted, however, that this additional source does 
not need to be a cluster group at the same redshift as 
C10016+16. It might be a chance alignment of two point 
sources, or it might be a group of galaxies at another red- 
shift. A single point-source is pretty unlikely to be the case, 
as we can reject a point-like shape with a confidence limit of 
99%. This confidence limit is determined by comparing the 
surface brightness distribution of the source with the surface 
brightness distribution of an artificial point-source having 55 
photons (see also FigJ^). The surface brightness distribution 
is marginally affected by the (subtracted) model parameters 
and the pixel size of the image. 

To exclude the possibility of an artifact produced by the 
satellite observatory (e.g. wobbling), we compare the quasar 
in the North of the cluster with a point source, again us- 
ing the surface brightness distribution. We obtain very good 
agreement of the quasar with a point-source, again the arti- 
ficial point source having the same number of photons than 
the quasar itself, after subtracting the cluster emission (note: 
the Poisson errors are calculated for the total number of pho- 
tons per bin, not only for the residual photons). Inspecting 
the optical image we find at a distance of about half an 
arcmin a red star. This star is by far the brightest optical 
source in this region. It is a M4V dwarf (BR). We can cal- 
culate an upper limit for the X-ray luminosity of the star, 
using the magnitude of the luminosity in the red band of 
m r =16.84, and m„=18.5 in the visible (from the Cosmos 
data base). The saturation limit for M dwarfs is found to 
have a value of 10 -3 (Fleming et al. 1993). This value de- 
fines the upper limit of the energy radiated in X-rays versus 
the energy emitted in the optical for M dwarfs. The 55 ex- 
cess photons found in the hard ROSAT band exceed this 
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Figure 7. The surface brightness profile of the additional source 
in the West of the centre of C10016+16. The crosses mark the 
data, inclusive errors, the dots show the distribution of a point 
source having the same number of photons as the original residual 
source. 

limit by a factor of 8, excluding this star as the source of 
the X-ray emission. 

However, we cannot exclude Quasars or AGN's to be the 
source, probably, if the case, in a chance alignment. 
A strong indication for the existence of a group of galaxies 
in the West of the centre of the cluster can be found by com- 
parison with the dark matter map obtained by SEFE using 
the weak lensing approach by Kaiser & Squires (1993). 

7.3 Comparison with weak lensing 

Recently SEFE presented a map of the dark matter distri- 
bution measured by weak lensing. They also made a rough 
estimation of the mass inferred by X-rays, also using the 
ROSAT/PSPC data. 

First of all, as SEFE already discussed, the PA of the 
mass distribution the light and the X-ray emission agrees 
well. Comparing our mass results with the SEFE mass esti- 
mates, however, yields a discrepancy. Our result of the pro- 
jected mass is 4-7x 1O 14 M0 integrating out to a projected 
radius of 600 kpc, whereas SEFE derive a projected mass of 
7.3x10 14 Mq using X-rays and a result for the weak lensing 
mass of 8.5xl0 14 Mq (depending on the mean redshift of 
the lensed background galaxies) out to the same distance. 
This discrepancy can be explained by using a different pro- 
jection model for the X-ray surface brightness, which is more 
peaked in the centre than observed with the limited resolu- 
tion of the PSPC data. Also shifting the lensed background 
galaxies to redshifts at l<z<1.4 can in principle resolve the 
discrepancies. The fact that the background galaxies might 
be at a higher redshift than 1 is also supported by an ob- 
servation of Luppino & Kaiser (1996). They measured the 
weak shear around MS1054-03, a rich and very X-ray lumi- 
nous cluster at z=0.83 (Luppino & Gioia 1995), and found 
a mass-to-light ratio in agreement with more nearby clus- 
ters, when the lensed background galaxies lie at a redshift 
of about 1.5. 

Comparing the morphology of the weak lensing mass map 
with our residual images yields very good overall similari- 
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ties. The bimodality which can be seen in the weak lensing 
maps exhibits similarities as seen in Fig.^, and also less strik- 
ing in Fig.|H| The structure in the West found in the PSPC 
data coincides well with the western extension in the map 
of SEFE. 



7.4 The gas to total mass ratio in C10016+16 

Integrating over the entire cluster (outer boundary 3Mpc), 
the gas to total mass ratio of the cluster is 14-32%. This is 
calculated for the Id fit. We did not take any uncertainties 
concerning the hydrogen column density into account. For 
the overall temperature we took 8.4 keV. These results are 
typical for clusters and undeline the effect of the Baryonic 
Catastrophe (Briel et al. 1992; White et al. 1993). The ratio 
does not vary very much with radius, it stays constant within 
the error bars. 



7.5 The Eccentricity of C10016+16 and its 
consequences 

In X-rays (ROSAT/PSPC and HRI data) C10016+16 shows 
an eccentricity of e = 0.18. As HRI and PSPC show the same 
eccentricity, the larger PSF of the PSPC does not seem to 
affect this value. It is in good agreement with the optical 
distribution of the galaxies of e ga i = 0.21 ± 0.02 derived by 
SEFE. The eccentricity is also comparable with results for 
other clusters like A2199 (Gerbal et al. 1992) and the sam- 
ple of Buote & Canizares (1992) (hereafter BC). BC's results 
on e are generally smaller than the value for C10016+16, but 
this might be caused by the limited area they took into ac- 
count for the determination (in all cases < IMpc). 
As we proof, the ellipticity of C10016+16 hardly affects the 
total mass profile of the cluster. The fact that the galaxy 
distribution and the X-ray surface brightness lead to similar 
eccentricities is a strong argument that we do not underes- 
timate the eccentricity of the dark matter, despite the fact 
that SEFE obtain a result for the eccentricity of the dark 
matter of e w i = 0.61 ± 0.06 (r « r c ) via the weak lens- 
ing approach of Kaiser & Squires (1993). This discrepancy 
might be caused by the limited area of the CCD's for the 
lensing approach. Therefore it might be, that SEFE only 
measure the eccentricity in the centre, which is likely to be 
higher than the rest of the cluster, due to the chain of galax- 
ies, and the results of BC, who found that clusters in X-rays 
tend to show a rounder appearance with increasing radius. 
It is also possible, that the eccentricity measured by SEFE 
is partly overestimated due to the clear bimodal shape of 
the dark mass map. Generally BC found out that X-rays 
are a good tracer for the shape of the dark matter, better 
than for example the galaxy distribution. However, the influ- 
ences on the gas mass profile are relatively high (deviations 
of up to 20%). Thus taking for the gas mass estimate only 
sectors of the whole cluster can bias the results. This is an 
approach X-ray astronomers sometimes undertake, to elim- 
inate regions which are strongly affected by substructure. 
However the differences are not drastic enough to change 
the gas-to-total-mass of the clusters dramatically. 



8 SUMMARY 

In this paper we present X-ray properties on C10016+16, 
a distant cluster of galaxies with a high X-ray luminos- 
ity of Ljf(2-10keV)=2.3-4.4xl0 45 erg/sec. The cluster has 
a temperature between 6-10keV. ASCA and ROSAT tem- 
perature results being well consistent. We present the total 
mass profile of the cluster based on X-ray measurements 
on ASCA and ROSAT, using a Monte-Carlo technique. The 
mass lies between 1.4-3.3 xlO 15 M@ integrated out to a ra- 
dius of 3Mpc. The gas-to-total-mass ratio lies between 14- 
32%, making this cluster another example for the so-called 
Baryonic Catastrophe (Briel et al. 1992; White et al. 1993). 
Applying a substructure analysis, we find indications for a 
not very massive group of galaxies falling into the centre of 
the cluster. This analysis is based on subtracting an ellipti- 
cal model cluster following the isothermal /3-model from the 
cluster X-ray image, with subsequent testing of the residu- 
als. This indicates, that C10016+16 is not very disturbed. 
The PA's of the distribution of the E and E+A galaxies, the 
weak lensing mass distribution, and the X-rays are coincid- 
ing very well with about 50°. 
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AG(x,y) = 



2na 



(SjJe-^'^-^'^m) 2 ^'dy') 1/2 (A2) 

To realize this function one can apply a modified Gauss 
filter on the original image with a filter size of 

and calculating the square root of this modified Gauss fil- 
tered image for each pixel and multiplying it with a factor 



of 2 ^ a ■ This gives AG(x,y). For the residual image the 
error AG(x, y) is the same as for the original smoothed 
image. Therefore to get the significance of the residuals 
one first calculates the subtracted count rate image by 
G(x,y)—Mode\\(x,y) and calculates the statistical error of 
the residuals which is AG(x,y). The significance image is 
then 

G{x,y) - Modell(a;, y) 



S(x,y) 



(x,y) 



(A3) 



APPENDIX A: THE SIGNIFICANCE OF 
SOURCES WITH A GAUSS FILTER 

In this appendix we describe how we determine the signif- 
icance of substructural features in a Gauss filtered image. 
The process of Gauss filtering can be described mathemati- 
cally by: 




f(x',y')dx'dy' (Al) 



x and y are the coordinates of the filtered and x' and y' 
the coordinates of the original image. The value of each 
pixel is described by f(x',y'). Applying the error propaga- 
tion formalism to this function, assuming Poisson noise (i.e. 
Af(x,y) = y/f(x,y)), we obtain: 
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